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ABSTRACT 

The paper deals with the one possible mechanism of the pulsar radio emission, i.e., with 
the collective curvature radiation of the relativistic particle stream moving along the 
curved magnetospheric magnetic field lines. It is shown that the electromagnetic wave 
containing one cylindrical harmonic exp{is(f>} can not be radiated by the curvature 
radiation mechanism, that corresponds to radiation of a charged particle moving along 
curved magnetic field lines. The point is that the particle in vacuum radiates the triplex 
of harmonics (s, s ± 1), so for the collective curvature radiation the wave polarization 
is very important and cannot be fixed a priori. For this reason the polarization of 
real unstable waves must be determined directly from the solution of wave equations 
for the media. Its electromagnetic properties should be described by the dielectric 
permittivity tensor e(w,k, r), that contains the information on the reaction on all 
possible types of radiation. 
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1 INTRODUCTION 

The curvature radiation is the type of the bremsstrahlung 
radiation when a radiated charged particle moves along the 
curved trajectory with the curvature radius po and its accel- 
eration is orthogonal to the velocity v. The cyclotron rota- 
tion of a charged particle in the external magnetic field B is 
the example of this motion when po = v±/uj c . Here uj c is the 
cyclotron frequency, uj c = eB/m e cy, e and m c are the charge 
and the mass of a particle, 7 is the particle Lorentz factor, 
and v± is the component of the particle velocity which is 
orthogonal to the magnetic field. 

Moving along the circular trajectory, a particle radiates 
at harmonics of the cyclotron frequency: oj = nuj c . This ra- 
diation is called as cyclotron radiation for a nonrelativistic 
particle and as synchrotron radiation for a relativistic parti- 
cle (7 3> 1). For the synchrotron radiation the maximum of 
the radiated power turns to be at large numbers of cyclotron 
harmonics: n ~ 7 s . The total radiation power / also grows 
with the particle energy, / oc 7 2 . Therefore, the synchrotron 
radiation of relativistic particles is presented widely in the 
space radiation (Ginzburg & Syrovatskii 1964). 

It is necessary to stress that the length of formation 
of the curvature radiation, though is larger than the wave 
length A, is much less than curvature radius po. So, the prop- 
erties of the curvature radiation do not differ from that of 
the synchrotron radiation in which the cyclotron radius is 
equal to the local curvature radius po. The frequency of the 
maximum of the spectral power, u) ~ cj 3 /po, and radia- 
tion power, / = 4/3e 2 C7 4 /poi increases with the particle 



energy. Here the dependence on 7 is stronger than for the 
synchrotron radiation because the curvature is fixed, and 
does not fall with the energy as for the motion in the con- 
stant magnetic field. 

The curvature mechanism of radiation is believed to be 
connected with the mechanism of the coherent pulsar ra- 
dio emission. Indeed, in the region of the open magnetic 
field lines in the pulsar magnetosphere there is a relativistic 
electron-positron plasma moving with relativistic velocities 
along curved magnetic field lines. For the typical values of 
curvature radius, po — 10 cm, and Lorentz factor of elec- 
trons and positrons, 7 ~ 10 2 , the characteristic frequency 
of the curvature radiation is in the radio band. In the mag- 
netosphere, where the curvature frequency coincides with 
the plasma frequency, cjp/7 3 / 2 , we can expect the collective 
curvature radiation. Here uj p = {4ne 2 n e /m e ) 1 ^ 2 is the usual 
plasma frequency, and n c is the plasma number density. In 
the strong magnetic field of the pulsar magnetosphere, when 
the charged particles can move along magnetic field lines 
only, the frequency of the plasma oscillations is 7~ 3 / 2 times 
less than the usual plasma frequency u p . 

It seems natural to continue the analogy between the 
curvature radiation and the cyclotron radiation for the col- 
lective radiation. But there is the essential difference be- 
tween them, which does not permit to rewrite formulas of cy- 
clotron plasma radiation for the curvature radiation replac- 
ing the cyclotron radius by the curvature radius. The matter 
is that at each point of plasma in the magnetic field the dis- 
tribution of particles over transverse velocities is isotropic. 
All directions of particle transverse motion exist, so that the 
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average velocity equals zero. It is not so for the curvature ra- 
diation when all particles have only one direction of motion 
along the magnetic field. 

For the plasma physics the problem of the collective 
curvature radiation is rather complicated since it demands 
the consideration of an essentially nonuniform plasma. It 
does not result from the change of parameters of magnetic 
field and plasma in space. These effects can be taken into 
account in the local approximation because the wave length 
of radiation is much less than the scales of inhomogeneities. 
In order not to lose the curvature radiation we need to in- 
clude into consideration the turn of the vector of anisotropy 
of the particle distribution function, /(p) oc 8(p — pyB/B), 
in space. Here py is the longitudinal particle momentum. 

Two parameters of the curvature radiation, i.e., the 
length of formation U — P0/7 — ^7 2 3> A and the width of 
the radiation directivity 8(f) — 7 _1 , connect by the relation 
If I po — 8(f>. Thus, the particle is in the synchronism with the 
wave (i.e., the particle sees the constant wave phase) along 
the path on which the wave intensity changes essentially. 
The value of A is the wave length of the curvature radiation, 
A ~ po/7 3 . 

The problem of calculation of the dielectric permittivity 
in the geometrical optics approximation for the nonuniform 
anisotropic plasma particle distribution was solved by Be- 
skin, Gurcvich and Istomin (below BGI, 1993). They also 
described the collective curvature-plasma interaction, when 
the electromagnetic waves, connected with the curvature 
radiation, are amplified simultaneously by the Cherenkov 
mechanism. This effect is absent in the vacuum. However, 
this procedure is rather complicated and demands clear un- 
derstanding. Because of that there are some incorrect state- 
ments in the literature (see, e.g., Nambu 1989; Machabeli 
1991, 1995). 

Apart, another way of investigation of the problem of 
the collective curvature radiation was carried out during 
many years (Asseo et al. 1983; Larroche & Pellat 1987; Lyu- 
tikov et al. 1999; Kaganovich & Lyubarsky 2010). They con- 
sidered more simple task connected with the pure cylindrical 
geometry which can be solved "exactly" . In such a statement 
the magnetic field lines are considered to be concentric, the 
relativistic plasma moving (i.e., rotating) along the mag- 
netic field lines owing to the centrifugal drift directed par- 
allel to the cylindrical axis (z-coordinate) with the velocity 
u — epe/po <S c. Here again p c = c/w c - But this approach 
cannot be used when analyse the curvature radiation (Be- 
skin, Gurevich & Istomin 1988). 

Indeed, let us choose the electromagnetic fields of the 
wave, as was done in all the papers mentioned above, in the 
form 

(E, B) = (E(p), B(p)) x exp {-iut + is(f> + ik z z} . (1) 

Here w is the wave frequency, s is integer number defining 
the azimuthal wave vector k^, and k z is the longitudinal 
wave vector along the cylinder. In this approach the wave 
amplitudes E(p),B(p) are to be considered as functions of 
the radial distance p only. Moreover, not vectors E and 
B, but their cylindrical components (E, B) p , (E, B)^ and 
(E, B) z depend on the coordinate p only. It means that the 
wave polarization follows the magnetic field, turning from 
one point (f> to another. It can be so if we have the defi- 
nite boundary condition, e.g., putting the system into the 



metallic coat. Under such suggestions we come to the one 
dimensional problem, which can be easily solved. Here we 
will show that such a wave does not have any relation to the 
curvature radiation. 

Really, let us consider the particle moving exactly along 
the circle of radius po with the constant velocity v; this mo- 
tion corresponds to the infinite magnetic field. Then the ra- 
diated power is equal to the work of the wave electric field 
under the particle electric current. The electric current is 



j = ev8((f> - tit) 8(z) 



<5(p-po) 



where f2 = v/po, and for selected polarization we get 



/ 



jEdr = evE t f > (po) exp {—iut + isflt} . 



(2) 



(3) 



As we see, the radiation is possible only if u — sQ, — 0, i.e., 
u) = k<pv. It is just the condition of Cherenkov, not curvature 
radiation. The point is that the wave with such polarization 
can not be radiated by the curvature mechanism. The differ- 
ence between the curvature wave and the Cherenkov wave 
is in the finite interaction time of the bremsstrahlung radi- 
ation with a radiated particle. The freely propagating wave 
with almost constant polarization deflects from the direc- 
tion of a particle motion. As a result, the nonzero projec- 
tion of the wave electric field on the particle velocity (i.e., 
on the direction of the electric current) occurs, and the wave 
takes away the energy from the particle. This continues the 
finite time r = lf/v that can be determined from the re- 
lation t(uj — kv) ~ 1. For the relativistic particle (v ~ c) 
r = (po/cjc 2 ) 1 / 3 ~ po/cy. Below we will find the real polar- 
ization of the curvature radiation. 

The paper is organized as follows. In section 2 we will 
find that the polarization of the curvature wave does not 
correspond to one cylindrical harmonic. In section 3 it is 
shown that the nonlinear wave interaction can lead to sig- 
nificant changes in cylindrical modes propagation. In section 
4 the BGI permittivity tensor will be derived from the per- 
mittivity corresponding to one cylindrical mode. Finally, in 
section 5 we discuss the main results of our consideration. 



2 POLARIZATION OF THE CURVATURE 
WAVE 

The radiation field of the electric current density j and the 
electric charge density p e of the moving particle with the 
charge e is described by the retarded potentials (Landau & 
Lifshits 1975): 



■m. 

R 



<1> 



J Pe(t') 



dr. 



(4) 
(5) 



Here t' = t — R/c is the retarded time, R is the distance from 
the charge location at the time t' to the observer which has 
the cylindrical coordinates (p, (f>, z), 



R = [p 2 + z 2 + po - 2pp cos((f>' - (f>)] 
d> = tit'. 



1/2 



(6) 
(7) 
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After the Fourier transformation of potentials (|4|-([5| over 
the time we obtain 



1 

2^ 
1 

2^ 



A(t) exp{iojt}dt, 
exp{iojt}dt. 



(8) 
(9) 



It is convenient now to replace the integration over time t by 
the integration over the retarded time t' and then over the 
angle ((f)' — (f)). As a result one can obtain for the cartesian 
components (a;, y, z) of the vector potential A and the scalar 
potential $ 

[A w ; *„] = |^ exp{iujcf>/n} \-K 3 ; K c ; 0; -7v l . (10) 

Z7TC L « J 

Here the quantities Ko, K s , K c are the functions of coordi- 
nates p and z only and they are equal to 



Ko = 



Ks 



da, 



K c 
R 



exp{iu(R/c + f2 _1 )a} 

R + vpsma/c 
cxp{iu)(R/ c + f2 _1 )a} sina 

R + vp sin a/c 
exp{iaj(i?/c + fl~ 1 )a} cos a 



R + up sin a/c 
(p 2 + z 2 + pi — 2pp cos a 



da, 



da, 



(11) 



,1/2 



The expression (10 1 is valid at any point r, i.e., not only in 
the wave zone. The dependence over the angle <f> is given by 
the exponent exp{iuxf>/Q,} . From the periodicity over (f> we 
have ui = sQ. 



The key point of the above expansion ( 10 \ is that the 
radiated wave is the superposition of three harmonics: s,s— 1, 
and s + 1. For example, the azimuthal electric field E^ is 
equal to 



Em 



lid 

V 



P c 



2 e 



Pa 



Kn 



(K s sin + Kc cos < 



.(12) 



The first term in Eqn. (12), which is proportional to the 



scalar potential $, is not important in the wave zone, 
p > po, but is significant in the near zone on the parti- 
cle trajectory p = po. Due to this term, the particle, which 
is in the resonance with one of three harmonics, say with 
s (u — sQ), is beaten out of the synchronism by neigh- 
bour harmonics s ± 1. The electric field E^ changes its 
sign during the time r. The synchronism condition, i.e., 



1 — cos Or : 



2/2 

V /c 



; 7 



defines the time r, 



r ~ 1/Q-y = po/cj, 



(13) 



which coincides with the time of formation of the curvature 
radiation. 

Thus, the radiated curvature wave consists of three har- 
monics s, s ± 1 with the fixed relation between their am- 
plitudes. Namely, this circumstance provides the curvature 
mechanism of the radiation. Appearance of harmonics s ± 1 
except the resonant one s = u/fl is due to the additional 
modulation of the radiation field induced by a modulation 
of the particle electric current having the harmonic s = 1. 
Now one can understand why the simple problem of the col- 
lective curvature radiation in the cylindrical geometry with 
only one azimuthal harmonic exp{is0} does not reveal any 



significant amplification of waves (Asseo et al., 1983; Lyu- 
tikov et al., 1999; Kaganovich & Lyubarsky, 2010). In this 
case the chosen wave polarization does not contain primor- 
dially the curvature mechanism. 



3 COLLECTIVE TRIPLE RADIATION 

In the previous section it was shown that the curvature radi- 
ation of one charged particle can not be described in the pure 
cylindrical geometry by one azimuthal harmonic exp{is<j!>}. 
In a collective radiation the modulation of the particle elec- 
tric current appears together with electromagnetic field ex- 
citation. Because of that the resonant azimuthal harmonic 
s — ujp/vtt, mixes with harmonics of the electric current mod- 
ulation and produces all possible values of s. Further in the 
section 4 we will see the all azimuthal harmonics s give con- 
tribution to the response of a media on an electromagnetic 
field. But in this section it will be demonstrated that the 
collective curvature radiation of only triplex of azimuthal 
harmonics (s, s ± 1) differs significantly from that of one 
harmonic s as it is usually considered in the literature. 

Let us consider the simple cylindrical one-dimensional 
problem of radiation of the cold stream of plasma particles 
with the charge e and the mass m e moving along the infinite 
azimuthal magnetic field Bo = B^. In this case the particles 
can move only in (^-direction with the velocity v$ at different 
cylindrical radius p. The unperturbed particle density n^ ' 
and velocity are constants, i.e., they do not depend on 
p. The electric current j has only ^-component as well as 
B 2 -component of the wave magnetic field (B p = B,f, = 0). 
Accordingly, the wave electric field has two components E p 
and E4, (E z = 0). 

The dependence of the wave fields over time and coor- 
dinates is the following 



[E p ; E^; B z ) = \E p (p); E^p); B z (p)] exp{-iut + is<j>}. 



Then, we obtain from Maxwell equations 



jpW . 2 



E 



("I 



dp 



a 



dE, 



00 



p * u p J $ 



E, 



00 



(14) 



(15) 



(16) 



dp p i" u) p"* 1 p 

Here index a corresponds to one of three harmonics s or 
s ± 1. For simplicity we use here the dimensionless variable 
r defined as r — piv/c, as well as quantities A = cjp/(o; 2 7 3 ) 
and Jo- = 4nj^ /(Aoj). Here oj p = (47me 2 /m e ) 1//2 is plasma 
frequency, and 7 is the Lorentz-factor of the particle motion: 
7 = (1 — i^/c 2 ) -1 / 2 . After these definitions the equations 
above take the following form 



dE 



00 



dr 



^OO-^OO 
r a 



r(<0 



dE, 



00 



dr r v r 



E, 



00 



(17) 



(18) 



As was already stressed, we consider here the interac- 
tion of three waves s, s ± 1. It is important that they are 
not independent and their interaction is realized by the 
static electric field [E p (p)\ E^,(p)]exp{icj>} having the first 
azimuthal harmonic s — 1. This electrostatic field turns 
to be the result of nonlinear interactions of high frequency 
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neighbour harmonics s and s ± 1. Equations for the mode 
s = 1 under the same definitions are 

AEj, i E$ 



dr 
AE P 
dr 



r 



-i-E^+AZ- 
r 



E„ 



(19) 
(20) 



Here Z = 4nnec/(Au>). 

To determine the response of the stream on the electro- 
magnetic fields of the wave one can use the continuity and 
Euler equations 

On 



dt 



■ Vfnv) 



0, 



(21) 
(22) 



It is easily to understand that only the (^-component of 
Euler equation is needed, while the radial component just 
provides us the equilibrium configuration across the infinite 
magnetic field. We represent the plasma number density and 



the plasma velocity as the expansion over powers of the wave 
amplitude 



v 4 , = vf ) +Sv^+Svf ) + ..., 
n = n (0) +Sn (1) + 8n (2) + .... 

The linear response can be easy found 



(i) 



^(0) kv 4> 
uj — kv 



(0) ■ 

<t> 



eE 6 



m c 7 3 (o; — fct>i°') 



(23) 
(24) 

(25) 
(26) 



where k = s/p. On the other hand, for the nonlinear current 
the nonlinear relation between Sv$ and 5p^> should be taken 
into account 



r 3c 3 (0) 5 ( Sv <i>) 2 

Sp^ = m e 7 Sv^ - -m e v\ ,r y y—. 

2 v c z 



(27) 



The result of cumbersome but straightforward calculation is 



Js = 



l-Svf/T 



e: 



1-svf/r 



+ a "(oy I — ; — (o), 

v\> \ 1- (s-M'/r 



"<t> 



El +1 El 



Js-l 
Js + 1 



l-{s-l)vf/r 



l-{s + l)vf/r 
Ei +1 E: 



E 



s - 1 



l-{s-l)vf/r 



r 

To) ' 



i-(s+i)v { ;>/ 

ElE 



1 (0) / 

l-svy/r_ 



r . E^E^ 



1 I i 1 \ (0) / ' " (0) " ' aT± . (0) / 

i-(s + i)vy/r v y i-svy/r_ 



4. ^<t> 



+ 



E%E, 



( s -i)* 



l l/r \ , , _ , . _l. i ,„"" /,w i _ ' i i /,.„ i _ ,.. _ n„ fl11 



(1 - (« + i)<70(i - «#7»0 (i - <V0(i - (* - 1)<70 

1 1 



(28) 
(29) 
(30) 
(31) 



i ■ (0) I 1 ■ (0) I 

1 — ivl jr 1 — jv\ jr 
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Here a = e/(m e c-y 3 uj) is the particle velocity divided 
over the velocity of light. The same quantities for plane 
waves can be found in (BGI, 1993). Equations above are 
evaluated with vacuum initial condition for the normal mode 
that can be presented analytically, = — J CT (r), Er = 

iaJ a {r)/r. Here Jo-(r) is the Bessel function. It should 
be noted that the singularity in equations (17 1 is passed 



smoothly by additional small term +ie in the resonance de- 
nominators in (1 25b (|26|). 



In numerical calculations equations (17l-(20l for a = s 



and a = s ± 1 were solved with two different values for the 
quantities J CT and Z. In the first case we neglect non-linear 



terms in ( 28 ) — ( 31 1 , while the second one corresponds to the 



full non-linear problem. On Fig. [T]the results obtained for 
this cases are presented. For better representation of the in- 
fluence of the nonlinear current, we choose the amplitudes 
of s — 1 and s + 1 modes twenty times higher than the am- 
plitude of the s mode. In reality the s-mode interacts with 
the whole continuum of modes, so this model assumption is 
rather reasonable. Fig. [I] shows that in this case the inten- 
sity of the wave \E\ 2 is approximate 2.5 times larger than in 



the case when the nonlinear current is neglected. Hence, one 
can conclude that three wave interaction is rather effective. 

Thus, we have shown that the triplex of cylindrical har- 
monics, which corresponds better to the curvature mech- 
anism, is amplified more effective than the separated har- 
monic having the single value of the azimuthal number. In 
fact the real polarization of the collective curvature mode 
can be obtained only by calculating the permittivity tensor 
of the streaming plasma in the strong curved magnetic field. 
The solution of wave equations produces not only the disper- 
sive equation for normal waves, oj = w(k), but defines also 
their polarization. A priory it is unclear what polarization 
corresponds to unstable modes. 

At first sight, the problem considered above is essen- 
tially nonlinear and has no direct connection with the ques- 
tion of the linear wave amplification. We included nonlin- 
earity only in order to connect harmonics s, s ± 1 self- 
consistently. And appearance of neighbour harmonics s ± 1, 
even for small nonlinearity, strongly change the s-mode am- 
plification. It is clear also that interaction of s ± 1 waves 
with the field s = 1 will result in all azimuthal harmonics. 
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Figure 1. Model calculations of two cases, A = 10 2 , v = 1GHz, 7 = 5, s = 125. 



4 TENSOR DERIVATION 

In this section we will show that the asymptotic behaviour 
of the BGI dielectric tensor in the case of large enough cur- 
vature radius po can be found directly from the plasma re- 
sponse on the one cylindrical mode. For the infinite toroidal 
magnetic field only the response to the toroidal component 
of the wave electric field E$ is to be included into considera- 
tion (Beskin 1999). Here and below we consider the station- 
ary medium only, so the time dependence can be chosen as 
exp{— iuit}. Making summation over all cylindrical modes, 
one can write down 



D*{p,4>) = E+{p,<l>)- ]T E l/> (p,s)K(p,s)exp{is<l>}, (32) 

s — — 00 



and the transition to cartesian coordinate system one can 
obtain: 



D x = E x + 



2 71 



p'dp'd(j)' 



p' 



£ E^p'^'Wp-p') 



S — — OO 



K(p, s) exp{is(4> — (/)')} sin 0, (35) 



Dy = Ey ~ 2^ 



p'dp'dtf 



£ E4p',cj>')d(p-p') 



K(p, s) exp{is((f) — <f>)} cos (f>. (36) 

We choose the local coordinate system with the t/-axis di- 
rected along the magnetic field and the :r-axis, that is or- 
thogonal to it. From the above equations one can obtain the 
permittivity kernel components 



where 



e yv (r,r') = l-^ S(p-p')K(p,s)x 

S — — oc 

exp{is((f) — </)')} cos (f) cos <f> ; 



u j 



v<f, df 



(0) 



dp,/,. (33) 



Here f^(p<j>) is unperturbed distribution function. Making 
the Fourier transformation 



E +iP,«) = ^ / E<,,(p,<t>')eM-™<t>'W (34) 



e yx (r,r') = — — S(p - p')K(p, s) 

s— — oc 

exp{is((j> — 4>')} cos cisin cj> ; 
11°° 

e xv (r,r') = — — S (p - p') K (P, s ) x 

s— — oo 

exp{is(cj> — 4>')} sintricos^', 



X 



(37) 



(38) 



(39) 



6 Ya. N. Istomin, A. A. Philippov, and V. S. Beskin 



i i 

ElI (r,r') = l-^ J2 S( P ~p')K(p,s)x 

s— — oo 

exp{is(4> — <f)')} sin 4> sin <f>' , 
that provides the material relationship 

A(r)= f e t] (v,r')E 3 (r')dr'. 



(40) 



(41) 



It should be noted that the operator presented above 
satisfies the needed symmetry condition 



eij(r,r',u) = e 3l (r', r, -w) 



(42) 



(it is provided by the condition K(r, s, ui) = K(r, —s, —ui))- 
As it is well-known (Kadomtsev 1965; Bornatici & Kravtsov 
2000), it is this symmetrical form of permittivity tensor that 
is to be used for the calculation of components of the per- 
mittivity tensor e,j (u), k, r) 



eij(u,k,r) -> r) 



£i7'(w,£,f7)exp{-jk£}dr. (43) 



Here rj = (r + r')/2, £ = r — r'. It is important that the 
above tensor only describes correctly wave-particle interac- 
tion in inhomogeneous media with slowly varying parame- 
ters (Bernstein & Friedland 1984). 

Substituting now the kernel components, one can find 



e xx (oj, k, 17) = 1 



271 



d£ exp{ — ik£} - 



l»?-€/2| 

]T 5(|»7 + €/2|H^-£/2|)^(|T7 + £/2|, S ) x 

s= — oo 

exp{is(0 — cj>')} sin cj> sin 4>' , 

If 1 

e xy (uj,k,r]) = — J dg exp{-ik£} x 

OO 

53 5(|r, + C/2| - |r7 - £/2|)#(|r| + €/2|,a) x 



(44) 



2tt 



exp{is(<j!> — 0')} sin cos 0', 

d£exp|— ikfj- — — x 

ll-£/2| 



(45) 



8(\r, + £/2\-\n-e/2\)K(]r, + t/2\,8) 



x 



s — — OO 



e„ tf (a;,k,?j) = 1- — y d£exp{-ik£} _ 



exp{is(<j!> — </>')} cos sin (f> , 
1 

x 



(46) 



*(l»7 + £/2| - |i| - £/2|)K(|i7 + £/2|, 8 ) x 

3— — OO 

exp{is(0 — (j)')} cos cos 0'. 



(47) 



In this equations, the angles cj> and <f>' are the functions of 
polar angles of vectors r) and £, a v , and 



cos 0' 



ry| sin a,, + (|£|/2) sina e 

cos a v — (|£|/2) cos 
" \V~U2\ ' 



(48) 
(49) 



As a result, integrals above are reduced to integration over 
£, which is perpendicular to rj. On the other hand, the ex- 



pression for the delta-functions in ( 44 1—| 47 1 is the following: 
5(0-7r/2) 



+ 



(|t, + €/2|-|f,-«/2|)J 

SjO + n/2) 
(\v + e/2\-\v-t/2\Ye' 



(50) 



where 9 is the angle between vectors r\ and £. So, the integra- 
tion over angles can be done easily. Finally, from the tran- 
sition rj — > r, one can obtain cos a,, — > cosa r = 1. Hence, 



according to (50 1 (k£) = where kn is the component 



of the wave vector parallel to external magnetic field. 

The property of the absence of k± is very important, it 
provides the same symmetry as it was in the case of homo- 
geneous medium: Ey (— cj, — k, — B, r) = Eji(u>, k, B, r) (Is- 
tomin 1994). This result differs from one obtained by Lyu- 
tikov at al. (1999). In this work the importance of transfor- 
mation ( 43 1 is neglected. 



Using finally the Taylor expansion over |£| and the re- 
duction of resonant denominator to delta-function, one can 
obtain: 

1 



oj\ri + £/2\/v$ - s 

^(l*7| 2 + l£| 2 /4) 1/2 



-/(...)* 



ds. 



As a result, one can write down 



.8ttV 



— £yx — 



F ( K ) ^5 d P<t» 



ui dp,/, 



F'{k) 



.87r 2 e 2 



1/3 2/3 

Po V 

2/3 1/3 

Po V df 



,1/3 



df 
dp$ 

■dpt. 



dpj,, 



Here 



-t-oo 

F(k) = - ^ exp{i/rf + « 3 /3}dt, 



2(u - fc||^) 2/3 



,1/3 2/3 



(51) 

(52) 
(53) 
(54) 

(55) 
(56) 



prime means the derivative, and po is the curvature radius 
of magnetic field. 

Due to high enough curvature radius of field lines in the 
pulsar magnetosphere, one can use the asymptotic behaviour 
of F(k) for k > 1 



FM 



i 

7TK 



2i 



After integration by parts, the final result is 



V 



"ii 

7 3 P0" 3 



(57) 



(58) 



Here by definition ui = u> — kv, and the brackets <> denote 
both the averaging over the particle distribution function 
fe+ e- {P<t>) an d the summation over the types of particles: 



<(•••)>= E f(-)fZ 



_(p<p)dP4>- 



(59) 
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We see that the tensor above is just the BGI tensor, that 
leads to instability of the so-called curvature plasma modes. 
In the limit po = oo this tensor, as expected, tends to the di- 
electric permittivity of a homogeneous plasma. The nonzero 
components e xy , e yx and Se xx = e xx — l in the tensor Ei 



for the finite curvature are due to nonlocal properties of the 
plasma response on the electromagnetic wave in curved mag- 
netic field. The parameter of nonlocality (v\\/u))/po is the 
ratio of the formation length of radiation to the curvature 
radius. For vacuum Gj ~ uj/j 2 , and the length Vn/ui coincides 
with the length of formation of the curvature radiation If. 

It is important that the components e xy = —e yx and 
Se xx essentially change the wave polarization. The relation 
between E$ and E p of the wave electric field, following from 
the tensor of the dielectric permittivity (|58[), is 



(e xy + n p n$) E4, + (8e xx + 1 - n%) E, 



0, 



(60) 



where n p and are components of the dimensionless wave 
vector: n = kc/w. For the tangent wave propagation (i.e., 
for n p = 0) we have E^, ~ (5e xx /e xy )E p ~ (c/ p c Cb)E p . As 
a result, a wave can produce the negative work under the 
electric particle current i.e., it can be excited. It is not 
so if 5e xx — e xy — when E<$, = 0. 



5 DISCUSSION 

Thus, as it was shown above, the wave polarization 
[E p (p); E#(p)] exp{is(j>} containing one cylindrical harmonic 
s suggests only the Cherenkov mechanism of radiation. In 
the curvature radiation mechanism of one particle in vac- 
uum the generated wave consists of three harmonics s, s ± 1. 
This property provides the exit from the phase synchronism 
of the wave with the particle motion which is inherent in the 
bremsstrahlung radiation. For the collective curvature radi- 
ation it is shown that the hydrodynamical model of plasma 
motion along the infinite magnetic field gives different re- 
sults of the wave amplification depending on the wave polar- 
ization. So, there is no another way to find the polarization 
of exiting waves instead of calculation of the response of the 
medium on an electromagnetic field, i.e., to use the dielectric 
permittivity. The correct procedure of dielectric permittiv- 
ity calculation using the expansion over cylindrical modes 
is also shown above. It was demonstrated that the tensor 
obtained by this procedure coincides with the BGI tensor 
calculated previously by another method. 

In conclusion it is worth to note that unsuccessful at- 
tempts to find the collective curvature radiation bring to the 



term 'curvature-drift instability' (Lyutikov et al. 1999). As 
was shown the chosen simple wave polarization, i.e., one s- 
harmonic, means only the existence of the Cherenkov mecha- 
nism of the wave generation. In this case the centrifugal par- 
ticle drift places the significant role. Practically all curvature 
effects come only to this drift. And the Cherenkov resonance 
on the drift motion produces small wave amplification in 
better case (Kaganovich & Lyubarsky 2010). Stronger mag- 
netic field produces less drift velocity and less Cherenkov 
effect though the curvature of a particle motion does not 
depend on the magnetic field strength at all. 
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